A full quantal theory of one-neutron halo breakup reactions 
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We present a theory of one-neutron halo breakup reactions within the framework 
^f) ' of post-form distorted wave Born approximation wherein pure Coulomb, pure nu- 

o ■ 

' clear and their interference terms are treated consistently in a single setup. This 



formalism is used to study the breakup of one-neutron halo nucleus ^^Be on several 
targets of different masses. We investigate the role played by the pure Coulomb, pure 
nuclear and the Coulomb-nuclear interference terms by calculating several reaction 
observables. The Coulomb-nuclear interference terms are found to be important for 
more exclusive observables. 
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I. INTRODUCTION 



■ Our current understanding of most of the processes governing the nuclear systems is 



based on studies made with stable nuclei, which constitute less than ten percent of all the 
nuclei known to exist in nature. Away from the valley of stability there are a large number of 
nuclei having very short half lives and very small one- and two-nucleon separation energies. 
Many of them exhibit a halo structure in their ground states in which loosely bound valence 
nucleon(s) has (have) a large spatial extension with respect to the respective core 0,0,0,0]. 
We still lack a fully microscopic understanding of the stability of these unique many body 
systems. 

These nuclei are also important from nuclear astrophysics point of view. The r-process 
and the s-process paths which together are dominant mechanisms for nucleosynthesis of 
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heavy elements above iron, pass through the region of neutron rich exotic nuclei, in the Segre 
chart. Properties of these nuclei are, therefore, important inputs to theoretical calculations 
of stellar burning which otherwise are often forced to rely on global assumptions about 
nuclear masses, decays and level structures extracted from stable nuclei. 



ing the structure and 
or 



Projectile breakup reactions have played a major role in unrave 
properties of halo nuclei. However it is clear that pure Coulomb js, |y, [?, 0, 0, [l^ 
pure nuclear [3, ll^ ll^ breakup calculations may not be fully sufficient to describe all the 
details of the halo bre akup data which have been increasing rapidly both in quality and 
quantity 12,Q, I3- In majority of them both Coulomb and nuclear breakup effects 



as well as their interference terms are likely to be significant and the neglect of the latter 



terms may not be justified 
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22|. A theory which can take care of the Coulomb 
and nuclear breakup effects as well as their interference terms on an equal footing is an 
important requirement in interpreting the data obtained from the experiments done already 
or are planned to be done in future. 

For breakup reactions of light stable nuclei, such a theory has been developed [2^ within 
the framework of post-form distorted wave Born approximation (DWBA), which successfully 
describes the corresponding data at low beam energies. However, since it uses the simplifying 
approximation of a zero-range 

is inapplicable to cases where the internal orbital angular momentum of the projectile is 
different from zero. 



Recently, we have presented a theory |26l . l27| to describe the breakup reactions of one- 
nucleon halo nuclei within the post-form DWBA (PFDWBA) framework, that consistently 
includes both Coulomb and nuclear interactions between the projectile fragments and the 
targets to all orders, but treats the fragment-fragment interaction in first order. The 
Coulomb and nuclear breakups along with their interference term are treated within a single 
setup in this theory. The breakup contributions from the entire continuum corresponding 
to all the multipoles and the relative orbital angular momenta between the valence nucleon 
and the core fragment are included in this theory where finite range effects are treated by 
a local momentum approximation (LMA) ■ Full ground state wavefunction of the 

projectile, of any angular momentum structure, enters as an input to this theory. 

The Coulomb-nuclear interference (CNI) terms have also been calculated using the prior- 
form DWBA and within models 0, 2^ where the time evolution of the projectile in 
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coordinate space is described by solving the time dependent Schrodinger equation, treating 
the projectile-target (both Coulomb and nuclear) interaction as a time dependent external 
perturbation. Recently, within an eikonal-like framework. Coulomb and nuclear processes 
have also been treated within the same framework in Ref. 

In this paper, we present more details of the formalism of PFDWBA breakup theory and 
of its application to breakup reactions of the one-neutron halo nucleus ^^Be on targets of 
masses spanning a wide range in the periodic table. We investigate the role played by the 
pure Coulomb, pure nuclear and the CNI terms by calculating different breakup observables. 
Our formalism is presented in section II. In section III, we present and discuss the results of 
our calculations for various observables for the breakup of ^^Be on various targets. Summary 
and conclusions of our work are presented in section IV. Additional discussions on the validity 
of the LMA relevant for our case are presented in Appendix A. 



II. FORMALISM 



We consider the elastic breakup reaction, a + t-^b + c + t, in which the projectile a 
[a = b + c) breaks up into fragments b and c (both of which can be charged) in the Coulomb 
and nuclear fields of a target t. The triple differential cross section for this reaction is given 
by 

,^jo^o = T^P^^»^ ^ l^^-l'' 
dEbdUbdUc hva ^ 

where Va is the relative velocity of the projectile with respect to the target, i is the orbital an- 
gular momentum for the relative motion of b and c in the ground state of a, and p{Eb, fi^, Qc) 
is the appropriate phase space factor (see, e.g., Ref. P]). The reduced transition amplitude, 
in Eq. ((T)), is defined as 

iPem(i^b,K;K) = J dvidvixt^*ikb,T)xir^*{K,rc)Vbciri) 

xuein)Yiih)xi^\K,r,), (2) 

with i = \^2£ + 1. In Eq. 0, functions Xi represent the distorted waves for the relative mo- 
tions of various particles in their respective channels with appropriate boundary conditions. 
The superscripts (+) and (— ) represents outgoing and ingoing wave boundary conditions, 
respectively. Arguments of these functions contain the corresponding Jacobi momenta and 
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coordinates. Vfec(ri) represents the interaction between b and c, and ue{ri) is the radial part 
of the c„.xesp„„di„g wavefu„c«o,U„ the ground state of a. The position vectors sat.sfy the 
relations (see also Fig. 1 of Ref . p ) : 

r = rj-ari, a = ■ , (3) 

rric + mi, 

rc = 7ri + 5r„ 6 = — ^ — , 7 = (1 - (4) 

where rrii [i = a, b, c, t) are the masses of various particles. 

The reduced amplitude [Eq- ©] involves a six-dimensional integral which makes its 
evaluation quite complicated. The problem gets further aggravated due to the fact that the 
integrand involves the product of three scattering waves that exhibit an oscillatory behavior 
asymptotically. In order to facilitate an easier computation of Eq. (j21), we perform a Taylor 
series expansion of the distorted waves of particles b and c about and write 

xl'\kt,r) = e-'^'^>'■^^x^~\Kr^), (5) 
x[r\K,r,) = e^^^-^^xi-Hkc,5r,). (6) 



Employing the LMA |28l . l29| , the magnitudes of momenta Kj are taken as 

K, (R) = ^{2m,/h\E,-V,{R)l (7) 

where rrij [j = b, c) is the reduced mass of the j — t system, Ej is the energy of particle 
j relative to the target in the center of mass (cm.) system, and Vj{R) is the potential 
between j and t at a distance R. Substituting Eqs. (0) and (0) in Eq. Q, the amplitude 
Pirni factorizes into two terms, each involving a three-dimensional integral. 



If X /, (8) 



where 



with 



and 



I rfrie-*Q-^^Hc(ri)ix,(ri)r4(ri), (9) 



Q = 7K, - aKb, (10) 



I = J rfraS+Vkb,r,)xi+)(-kc,^r,)xi+Hka,r,), (11) 
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where we have used the relation x^~^*(k, r) = k, r). 

Let us consider first the integral /. We expand the distorted wave for projectile-target 
relative motion in partial waves as 

xi^\K,r,) = T—i: ^'^fLSka,r,)Y^:*{K)YtSi^^), (12) 

where fiai^aj^i) is the radial part, calculated by solving the Schrodinger equation with 
proper optical potentials, which includes both Coulomb and nuclear terms. Beyond the 
range of the nuclear potential fiai^ai ^i) has the form 

hSka.rd "=°° ^e--«[4:^(A:ar.) - e^M^Kn)], (13) 

where H^j^'^ = Gl ± iFi, with and Gl being the regular and irregular Coulomb wave- 
functions, respectively, and is the scattering phase shift of the L^^ partial wave. In Eq. 
ai = argr(L + 1 + if]) is the Coulomb phase shift with Coulomb parameter rj. 
Similar expansions can be written for the distorted waves for the core fragment and the 
valence particle relative motions: 

X^Vk.,r.) = ^ E ^^''fLSh,n)Y^/,{k)Yt:{r.) (14) 

xi+)(-k„5r,) = E r^=^,(A;„5r,)F^,^(k,)F,^f (r,) (15) 

If ka(the incident beam direction) is chosen along the i-direction, then the spherical har- 
monic, Y^'^*(ka), in Eq. (fT^ simplifies to 

Y^f(^ = O,0 = O) = ^(5M„,o 

v47r 



^ Ma = O^Mb = -Mc = M (say). (16) 
Thereafter, substituting Eqs. (fT^. (fT^ and (fT3j) in Eq. (fTTj). we obtain 

/ = E {-r{^f-'^'-'^^UL.Y;i;/{h)Yt*{k) 

X {LbML,- M\LaO){L,OL,0\LaO) 

X / —fLa{K,ri)fLSh,ri)fLXkc,6ri). (17) 

Jo Ti 
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Let us now turn our attention to integral If [Eq. (jH))], which contains the structure 
information. Expanding the exponential, in Eq. (jU)), in partial waves, If simplifies to 

If = 4m-'Yi^{Q)J^ rldrMQri)u,{n)VUri). (18) 

Substituting Eqs. jlTj) and (HH)) in Eq. ©, we get 

X ytAh,k){LbOLcO\LaO)nL,,L.,LSka,h,ka), (19) 



where 



M 

oo 



Ze{Q) = / rfdnje{Qn)ue{n)VUri), (21) 
Jo 

nL,,L.,L. = / —fLAka,r.)fL,{h,r.)JUkc,Sr,). (22) 

The slowly converging radial integral IZl^^l^^l^ [Eq. (f^ ] can be effectively handled by 
using the complex plane method l^,!^- ^m^(Q)? Eq. (fTT?|l . whose argument contains 
the direction of a vector which is the sum of two other vectors [Eq. (jlUj) ]. can be expressed 
in terms of spherical harmonics corresponding to directions of those vectors as 3^ 

/ \ 1/2 

(IQI)X(Q) = E ^r^^ I i\c^K,\Y-'i^K,)' 

X {i - Lrrii - MLLMilimi) 

X F^;f^,jK,)r4(K,), (23) 

where L runs from to ^ and 

(24) 



is the binomial coefficient. 

This theory can be used to calculate breakup of both neutron and proton halo nuclei. 
Generally, the maximum value of the partial waves La, Li,, L^. must be very large in order to 
ensure the convergence of the partial wave summations in Eq. (|T9|) . However, for the case of 
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one-neutron halo nuclei, one can make use of the following method to include summations 
over infinite number of partial waves. We write f3em as 

r max 

Am = E MU) + E MU), (25) 

where (3 is defined in the same way as Eq. (jTHl) except for the summation sign and Lj 
corresponds to Lq, Lj,, and Lc. If the value of L™"^ is chosen to be appropriately large, the 
contribution of the nuclear field to the second term of Eq. ()25|1 can be neglected and we can 
write 

r max 

OO OO 

E A-(^«) ^ E /^fr'(^o - E /^fr'(^.), (26) 

where the first term on the right hand side, is the pure Coulomb breakup amplitude which 
for the case where one of the outgoing fragments is uncharged, can be expressed analytically 
in terms of the bremsstrahlung integral (see Ref. 0]). Therefore, only two terms, with 
reasonable upper limits, are required to be evaluated by the partial wave expansion in 
Eq. (ESI). 



III. CALCULATIONS ON ^^BE 

A. Structure model and optical potentials 

The wavefunction, U£{r), appearing in the structure term, Zi, has been calculated by 
adopting a single particle potential model in the same way as in Ref. The ground state 
of ^^Be was considered to be a predominantly s-state with a 2si/2 valence neutron coupled 
to the 0'^ ^°Be core [^°Be (g) 2s^/^2l'] with a one-neutron separation energy of 504 keV and 
a spectroscopic factor of 0.74 ^36|]- The single particle wavefunction was constructed by 
assuming the valence neutron-^°Be interaction to be of Woods-Saxon type whose depth was 
adjusted to reproduce the corresponding value of the binding energy with fixed values of the 
radius and diffuseness parameters (taken to be 1.15 fm and 0.5 fm, respectively). This gave 
a potential depth of 71.03 MeV, a root mean square (rms) radius for the valence neutron of 
6.7 fm, and a rms radius for ^^Be of 2.91 fm when the size of the ^°Be core was taken to be 
2.28 fm. The neutron-target optical pot entials used by us were extracted from t he g lobal 
set of Bechhetti-Greenlees (see, e.g, 3^), while those used for the ^°Be-target (jsj |3|) 
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system are shown in Table I. Following [31|, we have used the sum of these two potentials 
for the ^^Be-target channel. We found that values of L™""^' of 500 for Au, Ta, U, Pb and Ti 
targets and 150 for Be and C targets provided very good convergence of the corresponding 
partial wave expansion series [Eq. (fT!I|) ]. The local momentum wave vectors are evaluated 
at a distance, = 10 fm in all the cases, and their directions are taken to be same as that 
of asymptotic momenta (see appendix A). 



TABLE I: Optical potential parameters for the ^'^Be-target interaction. Radii are calculated with 
the rjt^/^ convention. 



system 


Vr 


rr 


ar 




ri 


ai 




(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) 


(fm) 


10Bg_197Au 


400 


2.08 


0.9 


76.2 


1.52 


0.38 


10Bg_208pb 


400 


2.08 


0.9 


76.2 


1.52 


0.38 


10Bg_44Ti 


70 


2.5 


0.5 


10.0 


1.5 


0.50 


i^Be-^Be 


100 


2.6 


0.5 


18.0 


2.6 


0.50 



B. Neutron energy distribution 

In Fig. 1, we present the results of our calculations for the double differential cross section 
as a function of neutron energy for two neutron angles (1° and 3.4°), in the breakup of ^^Be 
on Au at the beam energy of 41 MeV/nucleon. The core scattering angle in the laboratory 
system has been integrated from 0° to 30°. The dotted and dashed lines represent the pure 
Coulomb and pure nuclear contributions, respectively, while their coherent and incoherent 
sums are shown by the solid and dot-dashed lines, respectively. The plus signs and the 
inverted solid triangles represent the magnitudes of the positive and negative interference 
terms, respectively. 

The CNI terms are seen to be dependent on angles and energies of the outgoing neutron. 
Their magnitudes are nearly equal to those of the nuclear breakup contributions which leads 
to a difference in the incoherent and coherent sums of the Coulomb and nuclear contributions 
underlying thus the importance of these terms. 
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E„ [MeV] 



FIG. 1: Neutron energy distribution for the breakup reaction ^^Be on Au at the beam energy of 
41 MeV/nucleon, at the neutron angles of 1° and 3.4°. The dotted and dashed Unes represent the 
pure Coulomb and nuclear contributions, respectively while their coherent and incoherent sums 
are shown by the solid and dot-dashed lines, respectively. The plus signs and the inverted solid 
triangles represent the magnitudes of the positive and negative interference terms, respectively. 
The data are taken from [21 1. 

C. Relative energy spectra 

The relative energy spectrum of the fragments (neutron and ^°Be) emitted in the breakup 
of ^^Be on ^'^^Pb (top panel), ^^Ti (middle panel) and ^^C (bottom panel) targets at the beam 
energy of 72 MeV/nucleon is shown in Fig. 2. In these calculations the integration over the 
projectile cm. angle (6'„ioBe-Pb) has been done in the range of 0°-40°, mainly to include 
the effects of nuclear breakup coming from small impact parameters. The relative angle 
between the fragments (6'„_ioBe) has been integrated from 0° to 180°. The dotted and dashed 
lines represent the pure Coulomb and nuclear breakup contributions, respectively while their 
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FIG. 2: The differential cross section as a function of the relative energy of the fragments (neu- 
tron and ^''Be) in the breakup reaction of ^^Be on ^'^^Pb, ^^Ti and ^^C targets at 72 MeV/nucleon. 
The dotted and dashed lines represent the pure Coulomb and nuclear breakup contributions, re- 
spectively while their coherent and incoherent sums are shown by the solid and dot-dashed lines, 
respectively. The plus signs and the inverted triangles represent the magnitudes of the positive 
and negative interference terms, respectively. The data are taken from 

coherent and incoherent sums are shown by the solid and dot-dashed lines, respectively. The 
plus signs and the inverted triangles represent the magnitudes of the positive and negative 
interference terms, respectively. 

In case of breakup on a heavy target (^"^Pb) [Fig. 2 (top panel)] the pure Coulomb contri- 
butions dominate the cross sections around the peak value, while at larger relative energies 
the nuclear breakup is important. This is attributed to the different energy dependence of 
the two contributions jsil. The nuclear breakup occurs when the projectile and the target 
nuclei are close to each other. Its magnitude, which is determined mostly by the geometrical 
conditions, has a weak dependence on the relative energy of the outgoing fragments beyond 
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a certain minimum value. In contrast, the Coulomb breakup contribution has a long range 
and it shows a strong energy dependence. The number of virtual photons increases for small 
excitation energies and hence the cross sections rise sharply at low excitation energies. After 
a certain value of this energy the cross sections decrease due to setting in of the adiabatic 
cut-off. The coherent sum of the Coulomb and nuclear contributions provides a good over- 
all description of the experimental data. The nuclear and the CNI terms are necessary to 
explain the data at larger relative energies. 

In the middle panel of Fig. 2, we show the relative energy of the fragments in the 
breakup of ^^Be on a medium mass target (^^Ti). At low relative energies the pure Coulomb 
contributions are slightly higher than the pure nuclear ones, while at higher relative energies 
it is the nuclear part which dominates. Apart from the very low relative energy region 
the CNI terms play an important role, which is clearly borne out by the difference in the 
coherent (solid) and incoherent (dot-dashed) sums of the pure Coulomb and pure nuclear 
contributions. 

The relative energy spectra for the breakup on a light target (^^C) is shown in the 
bottom panel of Fig. 2. In this case we have used the same optical potential for the ^°Be- 
^^C system as in the ^°Be-^Be case, which we had used earlier in calculating the neutron 



angular distribution in Ref. 



26|. The total cross section in this case is normalized to the 



experimental cross section (found by integrating the area under the data points) and the 
same normalization constant is used for all the cross sections in this case. The breakup is 
clearly seen to be nuclear dominated at all relative energies, and the pure Coulomb and CNI 
terms have very little contributions. 

The importance of the peripheral region even in nuclear dominated reactions is underlined 
in Fig. 3. In this figure, we show a comparison of the angular distribution of the single 
neutron observed in the elastic breakup of ^^Be on a ^Be target at the beam energy of 41 
MeV/nucleon, calculated with (dashed line) and without (solid line) a lower cut-off of 10 fm 
in the integral [Eq. ()22|l ]. The close similarity of the two results shows the importance 
of the peripheral region in this reaction. This fact is further strengthened by noting that 
the total one-neutron removal cross sections calculated with and without cut-off are found 
to be 0.128 b and 0.189 b, respectively. Thus, almost 70% of this cross section is accounted 
for by regions larger than 10 fm. This is a consequence of the large spatial extent of ^^Be, 
which allows it to interact with the target nucleus even at a larger distance. 
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FIG. 3: The neutron angular distribution in the breakup of ^^Be on Be at 41 MeV/nucleon. The 
dashed hne shows the calculation where a lower cut-off of 10 fm is applied in the rj integral [Eq. 



H22|) ] while the solid line shows the result without any cut-off. The data are taken from 



D. Parallel momentum distribution 



The parallel momentum distributions (PMDs) of the ^°Be fragment in the breakup of 
^^Be on U and Ta targets, at 63 MeV/nucleon beam energy are presented in the rest frame 
of the projectile, in Fig. 4. The core transverse momentum has been integrated from 0-500 
MeV/c and the neutron angle (^„) has been integrated from 0° to 30°. We have used the 
following sets of optical potentials here: K- = 200 MeV, = 1.23 fm, ar = 0.9 fm, Wi = 
76.2 MeV, = 1.49 fm, = 0.38 fm for the U target and K = 200 MeV, = 1.26 fm, = 
0.9 fm, Wi = 76.2 MeV, = 1.53 fm, aj = 0.38 fm for the Ta target. The dotted and 
dashed lines show the contributions of the pure Coulomb and nuclear breakups, respectively, 
while their coherent sums are represented by solid lines. The coherent sum is normalized to 
the peak of the data, which are given in arbitrary units, and the same normalization factor 
has been used for the pure Coulomb and pure nuclear contributions. 

It is seen that around the peak region, the Coulomb contributions dominate. This is 
because most of the contributions in this region come from forward angles, where Coulomb 
breakup is the predominant mode. However, in the wings of the distribution (beyond about 
\pz\ = 40 MeV/c), contributions come from large scattering angles and consequently the 
pure nuclear breakup dominates, in this region. 
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FIG. 4: The parallel momentum distribution of the core in the breakup of Be on U and Ta 
targets, at 63 MeV/nucleon beam energy, in the rest frame of the projectile. The dotted and 
dashed lines represent the pure Coulomb and nuclear breakup contributions, respectively while 
their coherent sums are shown by solid lines. The data are taken from [s^. 

The FWHM of the distributions on U and Ta targets are found to be 48 MeV/c and 49 
MeV/c, respectively. The width of the fragment momentum distribution can be qualitatively 
related to the radial extent of the coordinate space nuclear wavefunction of the projectile 
via Heisenberg's uncertainty principle. Thus a narrow PMD width implies a large spatial 
extension of the nuclear wavefunction in the coordinate space for ^^Be. 

Our calculations with medium mass ^^Ti and light ^Be targets also leads us to similar 
conclusions. In Fig. 5 we present the PMDs of the ^°Be fragment in the breakup of ^^Be on 
^^Ti and ^Be targets, at 63 MeV/nucleon beam energy, in the rest frame of the projectile. 
The integrations over the transverse momentum and the neutron angle have been performed 
over the same range as in the case of heavy targets. The dotted and the dashed lines are the 
pure Coulomb and pure nuclear contributions, respectively, while their coherent sums are 
represented by solid lines. In case of breakup on a medium mass target (^^Ti) [Fig. 5 (lower 
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FIG. 5: The parallel momentum distribution of the core in the breakup of ^^Be on ^Be and ^^Ti 
targets, at 63 MeV/nucleon beam energy, in the rest frame of the projectile. The dotted and 
dashed lines represent the pure Coulomb and nuclear breakup contributions, respectively while 
their coherent sums are shown by solid lines. The pure Coulomb contribution in case of breakup 
on a ^ Be target is multiplied by a factor of 20 to make it visible. 

panel)] the pure Coulomb and pure nuclear contributions have nearly equal contributions, 
as has been seen while calculating the relative energy spectra of the fragments. Breakup 
on a light target (^Be) [Fig. 5 (upper panel)] is expectedly pure nuclear dominated. The 
pure Coulomb contribution is extremely small and is in fact multiplied by 20 to make it 
visible. The pure nuclear (dashed line) and the coherent sum (solid line) of the pure nuclear 
and pure Coulomb contributions thus almost coincide with each other. The FWHM of the 
distributions on ^^Ti and ^Be targets are found to be 60 MeV/c and 56 MeV/c, respectively. 
They are in line with the low widths mentioned earlier with heavy targets. 
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TABLE II: Total one-neutron removal cross section, various contributions from pure Coulomb and 
pure nuclear breakups, and their incoherent sum for ^^Be breakup on Au, Ti and Be targets, at 
beam energy of 41 MeV/nucleon. 



Target 


Total 


Pure Coulomb 


Pure nuclear 


Incoherent sum 


Expt. [21j 




(b) 


(b) 


(b) 


(b) 


(b) 


197Au 


2.10 


1.88 


0.21 


2.09 


2.5 ±0.5 


44 rp- 


0.403 


0.177 


0.189 


0.37 


0.55 ±0.11 


9Be 


0.189 


0.006 


0.181 


0.187 


0.24 ±0.05 



E. Total one-neutron removal cross section 

In Table II, we show the contributions of pure Coulomb and pure nuclear breakup mech- 
anisms to the total one-neutron removal cross sections in the breakup of ^^Be on Au, Ti 
and Be targets at the beam energy of 41 MeV/nucleon. The incoherent sum is obtained by 
simply adding the pure Coulomb and pure nuclear cross sections. 

For the heavy mass, high-Z target (Au) case it is seen that pure Coulomb breakup ac- 
counts for about 90% of the total cross section. On a medium mass, medium-Z target (Ti), 
the pure Coulomb and nuclear contributions to the total cross section, are nearly equal to 
each other, while for the low mass, low-Z target (Be), the pure nuclear contribution accounts 
for almost all of the total cross section. 

The total one-neutron removal cross section on Au and Be targets does not seem to be 
affected by the CNI terms, while for the Ti target case the incoherent sum seems to be 
about 10% less than the total cross section. Thus, it seems that the CNI terms manifests 
themselves more explicitly in more exclusive measurements, like double differential cross 
sections than in quantities like total cross sections. 

IV. SUMMARY AND CONCLUSIONS 



In this paper we have presented an extended version of a fully quantum mechanical theory 
of halo breakup reactions j^l within the framework of post-form DWBA, where the pure 
Coulomb, pure nuclear as well as their interference terms are treated consistently within 
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the same framework. In this theory, both the Coulomb and nuclear interactions between 
the projectile and the target nucleus are treated to all orders, but the fragment-fragment 
interaction is treated in the first order. The full ground state wavefunction of the projectile 
corresponding to any orbital angular momentum structure enters as an input to this theory. 
The lack of proper knowledge of appropriate optical potentials, particularly in the halo 
projectile-target channel is a source of uncertainty in the calculations. However this is the 
case for all reaction studies of halo nuclei where distorted waves in the projectile-target 
channel are used, as has already been pointed out in Ref. 3|. 

We applied our theory to study the breakup of one-neutron halo nucleus ^'^Be on sev- 
eral targets. Results for the neutron energy distribution in the breakup of ^^Be on Au 
at 41 MeV/nucleon emphasize the fact that the Coulomb-nuclear interference terms are 
both energy and angle dependent. They are almost of the same magnitude as the nuclear 
breakup contributions and this leads to a difference in the coherent and incoherent sums of 
the Coulomb and nuclear terms, more so at forward angles. The parallel momentum dis- 
tribution of the ^°Be fragment in the breakup reaction of ^^Be in the Coulomb and nuclear 
fields of U and Ta targets have also been calculated at 63 MeV/nucleon. It is seen that 
the region around the peak of the distribution, which gets substantial contributions from 
forward scattered fragments, is Coulomb dominated, while in the wings of the distribution, 
where contributions come from fragments scattered at large angles, the nuclear breakup 
contributions dominate. The FWHMs of the distributions were also found to be small con- 
sistent with the expectation of the wavefunction of ^^Be to have a large spatial extent in the 
coordinate space. Parallel momentum distributions in the breakup on medium mass (Ti) 
and light (Be) targets also confirmed a small FWHM. The relative energy spectra of the 
fragments (neutron and ^°Be) emitted in the breakup of ^^Be on Pb, Ti and C targets, at 
the beam energy of 72 MeV/nucleon have also been calculated. While the breakup on the 
light target was highly nuclear dominated, that on a heavy target required the nuclear and 
the CNI terms for a better explanation of the data particularly at higher relative energies. 
In case of breakup on a medium mass target, the total pure Coulomb and pure nuclear 
contributions were nearly equal in magnitude. The CNI terms were found to have little 
impact on the total one-neutron removal cross section in the breakup of ^^Be on heavy and 
light targets at the beam energy of 41 MeV/nucleon, but on a medium mass (Ti) target the 
CNI terms were almost 10% of the total one-neutron removal cross section. Thus in many 
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sophisticated experiments planned in the future one has to look into the role played by the 
CNI terms in analyzing the experimental data. 

The full quantal theory of one-neutron halo breakup reactions, presented in this paper, 
can be applied to describe the (a, 67) reaction provided the inelastic breakup mode is also 
calculated within this theory. These studies are in progress. There is also a need to extend 
the theory to describe the halo breakup at higher beam energies for which data have been 
taken at GSl, Darmstadt. This can be achieved by introducing the eikonal expansion of the 
distorted waves, instead of the partial wave expansion as done here. 
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APPENDIX A: VALIDITY OF THE LOCAL MOMENTUM APPROXIMATION 

The LMA provides a way of taking into account the finite range effects in the DWBA 
theory. It leads to the factorization of the breakup amplitu de [ Eg. Q], which makes its 
numerical calculation relatively simpler. As discussed in Ref. {28 1, a condition of validity of 
the LMA is that the quantity 

rj{r) = ^K{r)\dK{r)/dr\-^ (Al) 

evaluated at a representative distance R should be larger than the projectile radius (r^). In 
our case the LMA is done on both the outgoing fragments {b and c). In Fig. 6, we show the 
variation of rih{r) (top half) and r]c{r) (bottom half) with r, for the breakup of ^^Be on a 
Pb target, at the beam energy of 72 MeV/nucleon (left half) and ^^Be on a Be target at the 
beam energy of 41 MeV/nucleon (right half), where both Coulomb and nuclear potentials 
are included in the term Vj{r) in the definition of K{r) [Eq. ((7j)]. At r = 10 fm, r]b{r) = 185 
fm and rjclr) = 108 fm for ^^Be incident on a Pb target and ?7(,(r) = 914 fm and i]c{r) = 3892 
fm for ^^Be incident on a Be target. These values are much larger than the rms radius of 
^^Be, which is about 2.91 fm. (Incidentally, the spikes and turns in ri{r) at low r are due to 
the presence of the short ranged nuclear potential.) 
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FIG. 6: Variation of %(r) (top half) and r]c{r) (bottom half) with r, for the breakup of ^^Be on 
Pb, at 72 MeV/nucIeon (left half) and ^^Be on Be at 41 MeV/nucleon (right half). 

The variations of Kb[r) and Kc{r), the magnitudes of the local momentum (LM), with 
r for the above mentioned reactions are shown in Fig. 7. We see that Kh{r) and Kc{r) 
remains practically constant for r > 8 fm. Due to the peripheral nature of breakup reactions, 
this region contributes maximum to the cross section. Therefore, our choice of a constant 
magnitude for the local momentum evaluated at 10 fm is well justified. Thus the condition 
of validity of the LMA is well fulfilled in these cases. 

We have also performed calculations for different LM directions of b and c. We denote 
the different combination of directions as Di. both the LM angles of b and c are taken along 
asymptotic directions, D2: the LM angles of b are taken to be zero while the LM angles 
of c are taken along the asymptotic direction, D3: the LM angles of b are taken to be half 
those of the asymptotic direction while the LM angles of c are taken along the asymptotic 
direction, D^: the LM angles of b are taken along the asymptotic direction while the LM 
angles of c are taken to be zero, D^: the LM angles of b are taken along the asymptotic 
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FIG. 7: Variation of i^b(r) (top half) and Kc{r) (bottom half) with r, for the breakup of ^^Be on 
Pb, at 72 MeV/nucleon (left half) and ^^Be on Be at 41 MeV/nucleon (right half). 

direction while the LM angles of c are taken to be half those of the asymptotic direction, 
Dq-. both the LM angles of b and c are taken to be zero, and Dj: both the LM angles of 
b and c are taken to be half those of asymptotic directions. In Table III, we present the 

TABLE III: Calculated value of the total one-neutron removal cross section for ^^Be on Au and Be 
targets at 41 MeV/nucleon for different local momentum directions (see text). 

Projectile Di D2 D3 -D4 -D5 Dq D7 
+ target (b) (b) (b) (b) (b) (b) (b) 
^^Be + Au 2.10 2.13 2.14 2.16 2.10 2.24 2.29 
i^Be + Be 0.189 0.204 0.192 0.167 0.162 0.141 0.150 



total one-neutron removal cross section in the breakup of ^^Be on Au and Be targets at 41 
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MeV/nucleon for different directions of the LM. We note that these cross sections depend 
on the local momentum directions only to the extent of 10 — 15%. 




[deg] 



FIG. 8: The neutron angular distribution for the breakup reaction Be + Au Be + n + Au, 
at 41 MeV/nucleon for different LM directions (see text). 




[deg] 



FIG. 9: The neutron angular distribution for the breakup reaction Be + Be — >^ Be + n + Be, 
at 41 MeV/nucleon for different LM directions (see text). 

In Figs. 8 and 9, we show the variation of the neutron angular distribution for different 
combinations of LM directions of the core and the valence neutron in the breakup of ^^Be 
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on Au and Be targets, respectively, at 41 MeV/nucleon. In both the figures, the sohd fine 
shows the calculation with direction Di. In panel (a), of both Figs. 8 and 9, calculations 
are shown for directions D2 (dotted line) and (dashed line), panel (b) shows it for £^4 
(dotted line) and (dashed line), while panel (c) shows the same for Dq (dotted line) and 
Dy (dashed line). We see that for the Be target case the neutron angular distributions are 
minimally dependent on various LM directions, but for the Au target case, some dependence 
is observed at large angles. However, the magnitudes of the cross sections here are very small. 
Consequently no such angular dependence will be observed in cross sections which involve 
angular integrations for both the outgoing fragments or in more inclusive cross sections. 
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